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Comment on "Temporal scaling at Feigenbaum 
points and non-extensive thermodynamics" 

In a recent letter [1], P. Grassberger addresses the very 
interesting issue of tlie applicability of q-statistics to the 
renowned Feigenbaum attractor. However several points 
are not in line with our current knowledge, nor are the 
interpretations that he advances. 

To begin with, contrary to the statement in [1] there 
is a simple relation linking the constants in Eqs. (3) 
and (4). This is /3 = d/log2a, where is the fractal 
dimension of the attractor. The neglect of oscillations 
due to the multifractal nature of the attractor (his Ref. 
[31]) in the rate A„(l) leads to the above equality. The 
derivation of Eq. (3) for 2 > 1 is given in [2] so there is 
no need to reproduce it in [1]. 

A more important issue lies behind the Author's com- 
ment that Eq. (3) holds only for special values of time 
n. This is true as it is also true that there are many 
other special values of n that satisfy Eq. (3) exactly, all 
with the same value of q. See Ref. 36 in [1] and [3]. 
All together these sequences cover all n. Is there any key 
meaning behind this? As explained [3] , the dynamical or- 
ganization within the attractor is difficult to resolve from 
a simple time evolution: starting from an arbitrary posi- 
tion xq on the attractor and recorded at every n. What 
is observed are strong fluctuations that persist in time 
with a scrambled pattern structure. Conversely, unsys- 
tematic averages over xq and/or n would rub out the de- 
tails of the multiscale properties. However, if specific ini- 
tial positions with known location within the multifractal 
are chosen, and subsequent positions are observed only 
at pre-selected times, when the trajectories visit another 
region of choice, a well-defined g-exponential sensitivity 
appears, with q and the associated Lyapunov spectrum 
A(xo) determined by the attractor universal constants. 

Another point in case is the suggestion in [1] of adopt- 
ing Eq. (4) as focal point for the natural generalization 
of the Lyapunov exponent. This has already been con- 
sidered in Refs. [28] and [29] in [1], although, yet again, 
before taking a time average so that dynamical detail is 
preserved. A straightforward calculation shows that 


X{xo = 1) = 7^ In 
Inn 


dg 


dxQ 


- ln„ 


dg 


dxQ 


log2 a. 


where (1 — q)^^ — log2 a, and In^ y is the logarithm, 
the inverse of the g-exponential. So, the earlier definition 
for the generalized A is equivalent to that given for the 
same quantity by the g-statistics. The meaning of the 
index q is given by the above equalities. It is the degree 
of 'g-deformation' of the ordinary logarithm that makes 
A finite for large n. The physical origin of q is associated 
to the occurrence of dynamical phase transitions, of the 
so-called Mori's g-phase type, as demonstrated in [3]. 

The identity derived in Ref. [36] of [1] between the 
rate of g-cntropy change and the generalized Lyapunov 
exponent is not the identity {Sn — SQ)n~^ = A„n^^ in 
[1] (the zero identity for n — ^ 00) but refers to X{xo) 


as above. Of course it considers an instantaneous en- 
tropy rate K{xo) (comparable in the sense of [4] to the 
g-generalized KS entropy studied in Ref. [14] of [1]). The 
identity K{xo) — ^ixo) holds for n —^ 00 as the interval 
length (around xq) vanishes. It does fluctuate, but as 
explained, we look for the detailed dependence on both 
Xq and n. In contrast to the chaotic case there is not one 
identity but many, and the argument in [1] that aver- 
ages are needed for applications of Pesin's identity seems 
not to be useful for nonergodic and nonmixing trajecto- 
ries. Our results may not be insignificant as these can 
be reproduced [5] combining the arguments in Ref. [14] 
of [1] regarding the g-KS entropy with the results in [3]. 
One obtains the same equalities as for A(a:o = 1) above, 
with X{xq = 1) and \dg/dxo\ replaced by K{xo = 1) and 

C„(a;o)j respectively, where Cni^o) = 'zV^ and where 
Zn is Mori's partition function [3]. On the contrary, the 
Renyi entropies i/,^ in [1] from symbolic dynamics do not 
sense the universal constant a and/or the nonlinearity z. 

On the subject of the 'rich zoo' of g values, there is a 
well-defined family of these within the attractor, deter- 
mined by the discontinuities of the universal trajectory 
scaling function a [3] . There is a corresponding family of 
Mori's g-phase transitions, each associated to orbits that 
have common starting and finishing positions at specific 
locations of the attractor. The special values for g in 
the sensitivity are equal to those of the variable q in the 
formalism of Mori ei al at which the dynamical transi- 
tions take place [3]. Since the discontinuities' amplitudes 
diminish rapidly, there is a hierarchical structure in this 
family. The dominant discontinuity of a is associated to 
the most crowded and most sparse regions in the attrac- 
tor, and this alone provides a reasonable description of 
the dynamics for which the above expressions for A and 
K belong. About generality, a very similar picture has 
been recently obtained for another multifractal critical 
attractor, that of the quasiperiodic route to chaos [6]. 

A strong reason for preferring a g-exponential to a 
power law does not concern small arguments but the pres- 
ence of a time scale factor (the generalized A) absent (or 
hidden) in the power law. This useful quantity can be 
immediately 'read' from the anomalous sensitivity just 
like the ordinary A in chaotic dynamics. It is worth men- 
tioning that the (renormalization group) fixed-point map 
for intermittency, the other route to chaos, is rigorously 
given by a g-exponential map. (See Ref. [19] in [1]). 
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